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Abstract 

We perform an Amplitude Analysis of the world published data on 77 — * ir + n~ 
and 7r°7r°. These are dominated in statistics by the recently published results from 
Belle on the charged pion channel. Nevertheless, having only limited angular in- 
formation, a range of solutions remain possible. We present two solutions with 
r(/ (980) -> 77) = 0.42 and 0.10 keV, and r(/ 2 (1270) -> 77) = 3.14 ± 0.20 and 
3.82 ± 0.30 keV, respectively: the former being the solution favoured by x 2 , the 
latter at the edge of acceptability. Models of the structure of the /o(980) predict 
two photon widths to be between 0.2 and 0.6 keV, depending on its composition 
as mainly KK, ss or qqqq. Presently available data cannot yet distinguish unam- 
biguously between these predictions. However, we show how forthcoming results 
on 77 — > 7r°7r° can not only discriminate between, but also refine, these classes of 
partial wave solutions. 
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1 Introduction 



Two photon processes provide one of the cleanest probes of hadron structure and strong 
dynamics [lj. They are an almost unique way of measuring the charges carried by the 
constituents of any hadronic state. For neutral mesons built of quarks, a state R has two 
photon couplings given by 

r(a-ry) =1 E^n R | 2 , (l) 

q 

where the sum is over all the contributing quarks and in the non-relativistic quark model 
Ur is the probability that the constituent quarks annihilate. For heavy quarks this is 
proportional to the square of the wave-function at the origin. For light quarks there may 
be substantial corrections, estimated in Ref. [2] and reviewed in Ref. [3]. 
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Figure 1: Two photon decay rate of a meson in a quark picture is the modulus 
squared of the amplitude for 77 to produce a qq pair and for these to bind 
by strong coupling dynamics to form the hadron. 

Consequently, the two photon production of tt + tt~, for example, can be used to study 
the nature of all the states that couple to the system. Remembering that Bose sym- 
metry requires that odd J does not couple, we can access those states for which I and 
J are even. Among the most prominent of these is the /2(1270). Lying underneath this 
are the "enigmatic" scalars in the I = J = channel. Very many models have been 
proposed in the literature for their classification whether qq, qqqq or glueball — see for 
instance [3j HI [5] and the references therein. The two photon couplings of these scalars 
are one of the clearest guides to their compositions. However, the cr, / (980), and the 
many higher / ' s ! a H overlap with each other. This means that one is not looking to fit 
data with sums of Breit-Wigner forms and simply extract couplings to 77 of each reso- 
nance, but rather we have to perform a complete Amplitude Analyses to separate all the 
contributing spin components. Only then can we determine the couplings of the meson 
states by continuing these amplitudes to each resonance pole in turn. 

Data on the 77 — > 7r + 7r~ channel [6] have existed for many years, almost from the 
start of e + e~ colliders. Some of the earliest results, like that from DM2 [TJ, were merely 
event distributions, but with data from Mark II [8] at PEP twenty years ago followed 
by CELLO [10] at PETRA, normalised cross-sections for charged pion production 
became available. With the Crystal Ball detector first running at SLAC and then moving 
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to DESY, normalised cross-sections for 7r°7r° production in two photon collisions also 
resulted [H1[T2]- With these an Amplitude Analysis using a methodology that provides 
the nearest one can presently get to a model-independent partial wave separation of 
77 — ► 7T + 7r _ , 7r°7r cross-sections became feasible [HI [15] . This resulted in two distinct 
classes of solutions [16J, distinguished by whether the /q(980) appears as a "peak" or a 
"dip". Data in wide energy bins made both these possible. In the "peak" solutions the 
two photon width of the /o(980) was found to cover a range from 130 to 360 eV [H], as 
cited in the PDG tables [J?]. 

We now have data on charged pion production with a hundred times the statistics 
from Belle [TBI IT9"] allowing bins of 5 MeV width for both integrated and differential 
cross-sections. These reveal a clear peak in the 950-980 MeV region. Fitting the Belle 
integrated cross-section with model forms [19] gives a two photon width for the /o(980) 
of r(/ ) = 205lg3tJi7 eV. With the dramatically increased statistics of these results it 
is natural to consider how these improve the Amplitude Analysis, when the full data set 
(integrated and differential) are folded with the earlier world data. That is the aim of 
this work. 

We will find that the range of possible solutions is now more readily specified than 
in earlier work [T5l ITS] . With the large dataset from Belle, a two photon width for the 
/o(980) of 415 eV is favoured. However, solutions with a width from 100 to 540 eV are all 
acceptable. That such a wide range is still possible is partly because of the large systematic 
uncertainties in the Belle 7r + 7r~ angular distributions. However, these solutions can be 
resolved by precision 7r°7r° data particularly in key energy regions. To understand how we 
arrive at these conclusions, we outline the methodology in Sect. 2, present our solutions 
for 77 amplitudes solutions in Sect. 3 with plots of data fitted, and then in Sect. 4 discuss 
the results for resonance two photon widths, the implications for hadronic structure and 
give our conclusions. 



2 Methodology 

For a process, where one measures all observables and data cover the complete 47r angular 
range, such an Amplitude Analysis would be reasonably straightforward. However, here 
the photons are unpolarised and we have no information on the initial spin. Moreover, 
the detection of the charged final state in an environment of e + e~ means that pions 
are only identifiable in the centre of the angular range in the 77 rest system, typically 
I cos 9* I < 0.6. Consequently, one has large areas of the reaction plane with no direct 
experimental information. The effect of this is readily illustrated by an example. 

In Fig. 2 we show a compilation of the world data on 77 — > 7r + 7r~ from PEP, PETRA 
and BELLE [U US, IT9] . All are integrated over the range of | cos9*\ < 0.6. The dominant 
feature of the 77 — > nn cross-section below 1.5 GeV is the / 2 (1270) peak. This is a 
well-known spin-2 qq resonance with 1 = 0. It can be formed in two photon reactions 
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Figure 2: Comparison of the cross-section results for 77 — > tt + 7T~ from 
Mark II [8], Cello [9] and Belle [19]. In each case the cross-section is integrated 
over I cos#*| < 0.6. 

with either helicity or 2. The non-relativistic quark model [20] suggests pure helicity 2. 
Relativistic corrections have been found to be small [2]. What we observe is a large peak 
in the 77 — > mr cross-section. Let us assume we have picked out the total / = J = 2 
component. The angular distribution will then, in principle, tell us how much has helicity 
and how much helicity 2. In Fig. 3, the upper three plots show this distribution in 
a bin centred on 1270 MeV with different combinations of helicity and 2. If we had 
data covering the complete angular range up to |cos#*| = 1, the separation of these 
two helicity components would be relatively easy. However, if we now fold in the angular 
acceptance that cuts off at | cos 6 1 * | = 0.6, we arrive at the lower 3 plots. We see that 
these are virtually indistinguishable. To solve this problem, we need two things: firstly 
data of very high precision in the angular range covered. Secondly we need additional 
theoretical information to make up for the lack of complete angular coverage. Though the 
published results on 7r + 7r _ production from Belle [119] have sufficient statistical precision, 
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Figure 3: The upper plots show 3 examples of the angular distribution in the 
region of the /2(1270) with distinctly different ratios of I = S : D : D 2 
components. The lower plots show how these examples are typically modified 
when the angular efficiency fades away rapidly beyond cos#* = 0.6. 

systematic uncertainties mean that a large variation in the ratio of helicity 2 to helicity 
couplings of the / 2 (1270), and in the radiative width of the / (980), are still possible. 
As we shall see, precision 7r°7r° angular distributions, when finalised, are likely to resolve 
these ambiguities. With the present 7r + 7r~ data, the ambiguity displayed in Fig. 3 largely 
remains. The solutions with the larger two photon width for the /o(980) have a smaller 
Dq component, while the smaller S'-wave have a correspondingly larger D contribution. 
While the quark model for the / 2 (1270) would favour a small D contribution, and the 
Belle data [19] are consistent with this, our Amplitude Analysis does allow acceptable 
solutions with a larger D component. To see what progress can be made we need to 
understand the role of analyticity, crossing, unitarity and the low energy theorems of 
QED and chiral dynamics [21] . 

It is convenient to consider the Mandelstam plane for the process 77 — > tctt in the 
s-channel. This is shown in the upper plot of Fig. 4. The crossed-channels are, of course, 
77r — > 77r. At the threshold for this reaction, marked by the dot in Fig. 4, the photon has 
no energy and so has infinite wavelength. The probing photon then just sees the charge 
on the pions. Indeed, it is this Thomson limit of Compton scattering that defines the 
charge on the pion. At this point there is a low energy theorem due to Low [22J that 
requires the Compton amplitude to be equal to its Born term. This is given by one pion 
exchange in the t and ^-channels and a contact interaction to preserve gauge invariance. 
Thus at one point in this Mandelstam plane, we know the 77 — > 7r + 7r~ amplitude exactly. 
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Figure 4: Mandelstam plane with 77 — > 7r + 7r~ in the s-channel and Compton 
scattering off a charged pion in the t and w-channels. The (orange) dot marks 
the common threshold for the two Compton processes, where the Thomson 
limit applies. In the lower plot, drawn to scale, the shaded region of the s- 
channel delineates the physical region with —0.6 < cos#* < 0.6. The position 
of the t and w-channel pion poles indicating how very close they are to the s- 
channel physical region. The position of the p poles are shown, showing how 
far away the next most important crossed channel exchange are. The figures 
on the axes are in GeV 2 . 
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Abarbanel and Goldberger |23j showed this limit is approached in a smooth way so that 
the amplitude in the neighbourhood of the Thomson limit is still controlled by the one 
pion exchange amplitude. This is not really surprising since the pion mass is so very 
small. The pion poles at t and u ~ 0.02 GeV 2 are so much closer to the threshold than 
the contributions from any other exchange, such as that of the p with poles at t and u ~ 
0.5 GeV 2 . The position of these poles are drawn to scale in Fig. 4. In the lower plot we 
see the effect the experimental cut on cos 8* has on where we have data on 77 — > 7r + 7r~ in 
the s-channel physical region. The nearness of the pion pole to the forward and backward 
directions provides the key to how we can use the calculability of final state interactions, 
to make up for our lack of knowledge of the cross-section for cos 9* > 0.6 — at least at 
low energies, as we discuss briefly below. 

One important fact is that since the Born term that controls low energy 77 — > im has 
1 = 1 exchange in the t and u-channels, this involves both 1 = and 2 in the s-channel. 
Usually in reactions in which hadrons are produced in the final state, 1 = 2 amplitudes 
are generally much smaller than those with / = (or 1). This is because there are 
no established resonances with 1 = 2: indeed these are of course exotic in the quark 
model. Here in 77 — > titt at least at low energy, 1 = 2 amplitudes are equally important, 
otherwise the Thomson limit could not be maintained. This means that the reaction 
77 — > 7r + 7r - cannot be analysed in terms of amplitudes with definite isospin on its own, 
despite earlier attempts to do so. One must also have comparable information on the 
77 — > 7T 7r° channel too. Fortunately this comes from an earlier experiment by Crystal 
Ball at DESY [TT] with results limited to | cos 8*\ < 0.8. Our Amplitude Analysis treats 
both reactions simultaneously. The low energy theorem for 77 — > 7r°7r° is trivial, namely 
the amplitude is zero — neutral pions having no charge. Away from the Thomson limit, 
the chiral nature of pion interactions limits the scope for variation, as detailed in Refs. [T4"l 

HSH2HI3]. 

To proceed we first separate the amplitude for 77 — > tttt into its helicity components 
and then each of these into their partial waves. The unpolarized cross-section for dipion 
production by two real photons is given by the contribution of two helicity amplitudes 
M ++ and M H (the subscripts label the helicities of the incoming photons) : 

These two helicity amplitudes can be decomposed into partial waves as 

M ++ (s, 9, cf>) = e 2 VlQ^ £ F J0 (s) Yjo(6, cj>) , (3) 

J>0 

M + _( S , 9, 0) = e 2 v / 16V £ F J2 (s) Y J2 (6, 0) . (4) 

J>2 

The partial waves Fj\ (A = 0, 2) are the quantities we want to determine. 

Each partial wave Tj X is then an analytic function of s the square of the centre-of-mass 
energy E. Each has a right hand cut starting at the two pion threshold. Importantly, they 
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Figure 5: Cuts in the complex s-plane of the partial wave amplitudes for 
77 — > 7T7T, drawn to scale, illustrating how the crossed channel pion exchange 
controls the nearby part of the left hand cut. 

also have a left hand cut controlled by crossed-channel exchanges — Fig. 5. As we have 
emphasised, the part of this nearest to the direct channel is controlled by pion exchange 
as we have discussed. Corrections come from p, u, a%, ■ ■ ■ exchange, but these have a small 
effect on the 77 — ► tttt partial waves below 500 MeV. What is more important are tttt final 
state interations, especially in the I = J = channel. Dispersion relations [23] provide 
the vehicle for translating knowledge of this nearby left hand cut, the low energy theorems 
and experimental information on meson-meson scattering into precise information about 
the behaviour of the 77 — > tttt partial waves at low energy. This is described in Refs. [T4"l 
^ ne recen t work of Oiler et al. [26] confirms our calculation in the near threshold 
domain. A key aspect of this calculability of the final state interactions is that the 
discontinuity across the right hand cut is constrained by unitarity, as shown in Fig. 6. 
The sum on the right hand side is over all allowed intermediate states: allowed by quantum 
numbers and kinematics. At low energy only tttt contributes and the well-known Watson's 
theorem [29J applies. Since the 4tt channel is known to be weak, KK is the next most 
important contribution. Inclusion of just this channel allows our analysis to be extended 
up to 1400 MeV, then multi-meson states become critical as pp threshold is approached. 

If each partial wave for the reaction tttt — > hadronic final state H n with a given isospin 
is denoted by Tj , then coupled channel unitarity is solved by requiring for s > 4m^ 

?U S ) = «jaOO • T /( s ) > (5) 

where a(s) is an n-component vector of real coupling functions representing the coupling 
of 77 to channel H n . This reflects the fact that the right hand cut structure of the 
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Figure 6: Unitarity relation for each partial wave of 77 — > tttt. 

77 — > tttt amplitude is that of the hadronic processes with tttt final states. The functions 
a(s) have no right hand cutQ In general, there would also be a helicity label on the 
hadronic amplitudes too, but we will limit our consideration to spinless intermediate 
states like tttt and KK . Thus for the most part we will consider the energy region up to 
about 1400 MeV, for which a two channel representation is adequate, then more explicitly, 
we have 

^j,a(s; 77 -> tttt) = «ija(s) ?J{s) tttt -> tttt) + a 2 j X (s) fj(s; tttt -> TOsT) . (6) 

This is rendered useful by knowledge of these two meson-meson scattering amplitudes 
determined by combining experimental information, with unitarity and analyticity con- 
straints [32]. This provides a representation of the tttt — > tttt and 7r7r — > KK amplitudes 
from tttt threshold to 1400 MeV. The 'hat' indicates "reduced" amplitudes [33], in which 
for the S'-wave alone the pro cess- dependent Adler zeroes in the tttt — ► tttt and KK ampli- 
tudes have been divided outH Then Eq. (6) ensures that two body unitarity is respected 
for each 77 partial wave amplitude. The real functions a n (s) are the parameters in terms 
of which we fit the world two photon data. As already discussed, these are strongly con- 
strained (indeed even precisely known) down below 600 MeV by the Born term modified 
by calculable final state interactions. This serves to anchor our Amplitude Analysis at 
low energy. The use of the parametrization in Eq. (6) imposes continuity on our partial 
wave solutions. This helps to limit the ambiguity illustrated in Fig. 3 by relating what 
happens in one energy bin to that in its neighbours. In this way we perform an energy 
dependent amplitude analysis. With the latest data being in 5 MeV bins this is partic- 
ularly appropriate. This allows us to vary just 16-22 parameters to fit more than 2200 
data points as we discuss next. 



1 Very recently van Beveren and Rupp [30[ 131] have proposed that there is merit in writing these real 
functions in terms of particular complex functions. Though this can indeed be done, it introduces a right 
hand cut structure into these functions. The simplicity of our coupling functions a(s) is that they are 
real and having only left hand cuts are smooth for s > 4m^.. 

2 The I = J = hadronic amplitudes used [32] for the two-meson channels tttt and KK in the initial 
and final states have no zero in their determinant, and so the additional complications discussed in 
Refs. [331 [12] to parametrise the coupling functions is not required. 
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3 Amplitudes determined: solutions A and B 



Armed with the methodology of Sect. 2 and the representation of Eq. (6), we can fit the 
complete set of experimental cross-sections, integrated and differential, listed in Table 1 by 
expressing the functions Q\(s), a 2 (s), in terms of smooth functions exactly as in Refs. [15], 
[16]. In Table 1 we see that the number of data-points is far greater for the reaction 
77 — > 7r + 7r~ than for 77 — > tt°7[°. Since an accurate separation of the 1 = component 
requires both are accurately described, we give different weight factors to each dataset. We 
choose these so that the Mark II, CELLO and Belle data have roughly the same number 
of weighted data, while the weight assigned to the Crystal Ball data approximately equals 
the weighted sum of the 7r + 7r~ data. Nevertheless, good agreement is not easy to achieve. 

The analysis program works by integrating the amplitudes over the appropriate bin in 
energy and angle for each data-point. It does not just use the central values. This is to 
allow for any strong local variation of the amplitudes, particularly near KK threshold. In 
the many figures displaying the solutions we find, this should be borne in mind. Where 
the energy bins are sizeable (as with Crystal Ball) histograms are plotted (see Fig. 8). 
Where the energy bins are fine (as with Mark II data in 10 MeV steps and Belle data in 
5 MeV bins), the fits are shown more appropriately as continuous lines joining the bin 
centres, as in Figs. 7, 9, 10, 11 and 12. However throughout, the fits are histograms of 
appropriate bin width. 

Table 1: Number of data in each experiment below 1.44 GeV. Mark II results are 
from Boyer et al. [8], Crystal Ball from Marsiske et al. (CB88) [11] and Bienlein 
et al. (CB92) [12], CELLO from Harjes [9] and Behrend et al. [10], and Belle from 
Mori et al. [31]. 
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Figure 7: Solutions A and B compared with the Crystal Ball results on 77 — > 
7r°7r° . The 1988 Crystal Ball data [11] are integrated over |cos#*| < 0.8, 
while the 1992 data |12| with increased statistics cover I cos 0*1 < 0.7. 



The solutions are determined by the published statistical errors. In the case of the Belle 
results, the point-to-point systematic errors are folded into these statistical uncertainties. 
This still leaves the Belle results as the most precise of all the two photon data available, as 
illustrated in Fig. 2. However, only CELLO have folded their full systematic uncertainties 
into their published results. For all others there is a global systematic uncertainty that 
allows all the cross-sections across all energies to be renormalised. For Mark II this is 
by 11%, for Crystal Ball the shift is 11% for the 1988 results and by 3% for their 1992 
data. Only for Belle is the energy dependence of the systematics determined. This is 
typically 12%, completely dominating the statistical errors. We therefore quote the x 2 i n 
two ways for our solutions. We give that fixed by the "statistical" errors and in brackets 
that given by folding in the full global systematic uncertainties in quadrature. In our 
plots, displaying solutions against the data they fit, the overall systematic errors are not 
included. 
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The Mark II, Crystal Ball and Belle datasets [8j [TTJ [121 HH] each have a systematic 
normalization uncertainty. To fit the Amplitudes to the competing datasets, we allow each 
to shift within these bands. In the Solutions we display, rather than plot the published 
data shifted, we choose to rescale the fitted amplitudes. Consequently, the solutions 
displayed in Figs. 7-18 may appear slightly shifted between one dataset and another. 
This is most evident in Fig. 8, where the near threshold normalization of Mark II differs 
from that above 450 MeV. Though the Belle data [TU] have an overall 11-12% systematic 
uncertainty, the solutions generally take up no more than a 3% shift. 

In terms of x 2 , Solution A is favoured and provides a very good description of the 
world data. While the "statistical" errors give a y 2 per degree of freedom of 2.9 for this 
solution, with systematic uncertainties folded in this reduces to 0.9. This solution has 
a prominent /q(980) peak in the I — J — cross-section, and a smaller D component 
throughout the energy region covered. Consequently, this solution predicts a prominent 
/o(980) signal in the 77 — > ir°7t° cross-section. 

When systematic uncertainties are allowed for, Solution A has a probability of 99.9%. 
The range of other acceptable solutions is somewhat subjective, largely because of the 
complex systematics as well as inconsistencies between datasets, as seen in Fig. 2. The 
criterion we use is not precisely defined statistically, but based on the goodness of fit to the 
eye. We have explored the space by requiring varying amounts of D to D 2 components, 
as described in Refs. [TJl [TB]. Qualitatively, the solutions are either like Solution A, 
but with worse x 2 \ or closer to what we now describe as Solution B. This is the second 
solution we display. This solution has an overall x 2 of 1.2 per degree of freedom (being 3.7 
from wholly "statistical" errors). Though this solution only has a probability of 2.10~ 7 , 
it provides an acceptable fit to the Belle 7r + 7r~ data, when systematic uncertainties are 
folded in. It is the other datasets it mistreats. Solution B has a much smaller / (980) 
signal in the I = J = cross-section, while having a larger I = D component. The fit 
to the distinct peak in the Belle cross-section between 950 and 1000 MeV for 77 — > tt + tt~ 
in this solution is just as good as that in the favoured solution. The peak in the 7r + 7r~ 
cross-section, integrated over cos#* from —0.6 to +0.6, is reproduced by a constructive 
interference between the I — S and Dq components. 

Figs. 7-10 show how well Solutions A and B fit the integrated cross-sections, x 2 favours 
amplitude A, but B is an example at the edge of acceptability. In all cases the fitting 
of the Belle 7r + 7r~ results for the integrated cross-section between 800 and 900 MeV is 
not good — see Fig. 10 (and the comparison of data in Fig. 2). This is surely correlated 
with the strange cos 9* dependence in the Belle results [TH] shown in Fig. 15 in the same 
mass region. This is most likely generated by the strong sensitivity to the details of the 
separation of the large background in this experiment. 
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Figure 8: Solutions A and B compared with the Mark II results p3] on 77 — ► 
7r + 7r~ integrated over | cos#*| < 0.6. The freedom to make a small systematic 
shift in normalisation of the data above 0.4 GeV has been used to improve the 
fits. Rather than shift the data, the solutions have been renormalised. This 
results in the discontinuity in the solid curves at 0.45 GeV 



13 



Integrated cross-section yy — n + n 



300 



a (nb) 



200 



100 



i 1 r - 



A Cello 






0.2 0.4 0.6 0.8 1.0 1.2 1.4 



300 



a (nb) 



200 - 



100 







~ i 1— 



A Cello 




J i_ 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 

E (GeV) 

Figure 9: Solutions A and B compared with the CELLO results [9j [10] on 
77 — ► 7r + 7r~ integrated over | cos9*\ < 0.6. 
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Figure 10: Solutions A and B compared with the Belle results [19] on 77 
7r + 7r~ integrated over |cos#*| < 0.6. 
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Figure 11: Differential cross-section for 77 — > 7r°7r° from the Crystal Ball 
experiment [TU [12] compared with Solution A from the Amplitude Analysis. 
The numbers give the central energy in GeV of each angular distribution listed 
in order of the size of the cross-section at z — 0, where z = cos 9*. 



We now turn to the angular distributions. In Figs. 11-18 we compare Solution A 
with the differential cross-sections from Mark II, CELLO and Belle for charged pions 
and Crystal Ball (1988 and 1992 results) for tt°tt°. While the huge dataset from Belle 
dominates, there appears to be a structure in the | cos#*| ~ 0.45 — 0.55 region below 
1.1 GeV, seen in Fig. 15, to which we have just alluded. The fact that this possible 
"measurement bias" gives the fits a poor x 2 is one of the major reasons that the Amplitude 
Analysis presented here cannot be more precise about the /o(980) radiative width. This 
poorly fitted "structure" in the angular distributions makes a whole range of solutions 
equally acceptable. Solution B provides fits to the angular distributions that are not 
significantly different in the whole region below 1450 MeV, and so we do not show these 
fits separately. The somewhat worse x 2 apparent in Table 2 are barely discernible by eye. 

An important aspect of the Amplitude Analysis is that the individual partial waves 
are determined at low 77 energy by the QED low energy theorem and chiral constraints 
on the 7T7T final state interaction. These are implemented by the use of dispersion relations 
as discussed in Sect. 2 and as calculated in Refs. [23, EI]- The precision at low irn masses 
comes from the dominance of one pion exchange to cross-channel forces (as in Fig. 2), 
modified by calculable final state interactions. All solutions are required to respect these 
properties. As the energy increases other crossed-channel exchanges become increasingly 
important. As their contribution is less certain, the solutions are allowed to vary within 
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Table 2: Summary of contributions from each experiment to the total x 2 for solution 
A and B. Here x?ot * s calculated by dividing the sum of the x 2 ' s f° r ai l datasets by 
the total number of data-points we are fitting, namely 2267. XaVerage * s com Puted 
by dividing the sum of the x 2 's for each dataset by the total number of data-points 
in that experiment. * indicates that the x 2 includes statistical and point-to-point 
systematic errors. Only those with no asterisk include the overall systematic error. 



favoured SOLUTION A xL = 2 -88* (0.91) 


Experiment 


Process 


data-points 


V 2 

A average 


V 2 

A Int. X— sect. 


V 2 

A Ang.distrib. 


Mark II 


77 — > 7T + 7T~ 


144 


1.64* (1.02) 


1.52* (0.94) 


1.80* (1.12) 


Cr. Ball 


77 — > 7r°7r° 


126 


2.06* (1.72) 


2.04* (1.51) 


2.07* (1.81) 


CELLO 


77 — > 7T + 7T~ 


333 


1.76 


0.76 


1.57 

from Harjes 
2.10 

from Behrend 


Belle 


77 — > 7T + 7T~ 


1664 


3.34* (0.67) 


0.85* (0.12) 


3.55* (0.72) 




SOLUTION B xL = 3-72* (1.16) 


Experiment 


Process 


data-points 


A average 


v 2 

A Int. X— sect. 


v 2 

A Ang.distrib. 


Mark II 


77 — > 7T + 7T~ 


144 


2.01* (1.17) 


1.76* (1.06) 


2.32* (1.31) 


Cr. Ball 


77 — > 7T 7T 


126 


2.63* (2.27) 


1.95* (1.45) 


2.91* (2.60) 


CELLO 


77 — > + 


333 


2.10 


0.89 


1.94 

from Harjes 
2.22 

from Behrend 


Belle 


77 — > + 


1664 


4.33* (0.89) 


1.00* (0.19) 


4.61* (0.95) 
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Figure 12: Differential cross-section for 77 — > 7r + 7r~ from the Mark II ex- 
periment [8] compared with Solution A from the Amplitude Analysis. The 
numbers give the central energy in GeV of each angular distribution listed in 
order of the cross-section at z = 0, where z = cos 9*. The data are nor- 
malized so that the integrated cross-section is just a sum of the differential 
cross-sections in each angular bin. 

the band of this uncertainty. In Fig. 19 are shown the bands for the I = S, D , D 2 
waves and how Solutions A and B fit within these. The corresponding 1 = 2 waves are 
similarly controlled by the meson-exchange amplitudes, but as they have much weaker 
final state interactions, they are closer to their Born amplitude. These are not shown. 
All waves with J > 4 are approximated by their Born contribution in the whole energy 
range fitted, i.e. up to 1.45 GeV. They generally make up only a small component of the 
7r + 7r~ and 7T 7r° amplitudes and so this is a reasonable guide to their contribution in the 
energy range covered. 
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Figure 13: Differential cross-section for 77 — ► 7r + 7r~ from the CELLO ex- 
periment pj compared with Solution A from the Amplitude Analysis. The 
numbers give the central energy in GeV of each angular distribution listed in 
order of the cross-section at z — 0, where z = cos 9*. 
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Figure 14: Differential cross-section for 77 — ► 7r + 7r~ from the CELLO exper- 
iment in finer binning in both energy and angle from those of Fig. 11 [TO] , 
compared with Solution A from the Amplitude Analysis. The numbers give 
the central energyin GeV of each angular distribution listed in order of the size 
of the cross-section at z = 0, where z = cos 9*. 
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Differential cross-section 
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Figure 15: Differential cross-section for 77 — ► 7r + 7r~ from the Belle experi- 
ment [19J compared with Solution A from the Amplitude Analysis. Each graph 
in 25 MeV intervals displays data and the result of the Amplitude Analysis in 
5 MeV bins, so each graph has five sets of data and 5 curves, z = cos 9*. 
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Figure 16: Differential cross-section for 77 — ► 7r + 7r~ from the Belle experi- 
ment [19] compared with Solution A from the Amplitude Analysis. Each graph 
in 25 MeV intervals displays data and the result of the Amplitude Analysis in 
5 MeV bins, so each graph has five sets of data and 5 curves, z = cos 6*. 



21 




Figure 17: Differential cross-section for 77 — > 7r + 7r~ from the Belle experi- 
ment [19] compared with Solution A from the Amplitude Analysis. Each graph 
in 25 MeV intervals displays data and the result of the Amplitude Analysis in 
5 MeV bins, so each graph has five sets of data and 5 curves, z = cos 9*. 
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Figure 18: Differential cross-section for 77 — > 7r + 7r~ from the Belle experi- 
ment [19] compared with Solution A from the Amplitude Analysis. Each graph 
in 25 MeV intervals displays data and the result of the Amplitude Analysis in 
5 MeV bins, so each graph has five sets of data and 5 curves, z = cos 9*. 
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Figure 19: The bands display the moduli of the three dominant 1 = partial 
waves S, D and D 2 at low energies as determined by the dispersive analysis 
of the QED and chiral constraints discussed in the text. These amplitudes are 
the appropriate Fj x defined in Eqs. (1-3). Precision close to threshold comes 
from the dominance of one pion exchange to the cross-channel forces — see 
Fig. 2. For orientation the 1 = Born S'-wave (labelled B) is the black line. As 
the energy increases p, cu, a±, exchanges become increasingly important. 
Their contribution is less certain. This is reflected in the width of the bands. 
The solid lines are these same moduli in the Solutions A and B. 
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4 Two photon widths of the low mass isoscalars 



We are now in position to show the partial wave content of our Solutions A and B. In 
Fig. 20 we illustrate the contribution to the integrated cross-section of each of the 1 = 
and I = 2 S, D , D 2 waves and I = 2 S and D 2 waves. Very close to threshold these 
are controlled by the one pion exchange Born term. If this were not modified by different 
final state interactions in the 1 = and 1 = 2 channels the 1 = cross-sections would 
be exactly a factor of 2 larger than those with 1 = 2. One sees the difference is greater 
reflecting the stronger tttt interactions in the isoscalar channel. These bring about the 
marked dip in the 1 = S-wave seen around 550-600 MeV in Figs. 20, 21. The sizeable 




Figure 20: Contributions of the individual partial wave components to the full 
integrated cross-sections for Solutions A and B. Only the dominant waves are 
shown. The 1 = 2 amplitudes are presumed to be computable from the Born 
amplitude modified by known final state interactions, and so are the same for 
all solutions — as for the 1 = amplitude close to threshold. 
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difference between the Born amplitude and this S-wave is the signal of the cr-resonance, 
as noted in Ref. [25J. As the energy increases the isoscalar scalar component grows with 
a clear peak of the /o(980) in Solution A, and a little less so in Solution B. It then has a 
significant contribution above 1200 MeV that one would associate with the /q(1370). The 
largest component above 500 MeV is, of course, the / 2 (1270). In the favoured Solution A 
this appears overwhelmingly in the D 2 wave, while in Solution B this is shared with the 
D Q wave. To fit the data with limited angular information one has the obvious ambiguity 
of a significant S-wave in Solution A or a more complicated interference between S and 
Dq components in Solution B. Present data favour Solution A, but solutions like B are 
not yet ruled out. 

In Fig. 21 are shown the Argand plots for the important 1 = S-wave amplitude 
for Solutions A and B. This S'-wave amplitude encodes the effect in experiment of the 
a, /o(980) and /o(1370) poles, the first and third of which are very wide. Consequently 
the amplitudes vary greatly from the real axis, where experiment is performed, to the 
pole positions. The Argand amplitudes are seen to move fastest between 950 and 1025 
MeV, as expected from the far narrower / (980). Solution A shows a larger movement 
indicating the larger radiative width for the /o(980) in this solution. Fig. 21 also shows 
the slowly changing Argand plot for the 1 = 2 S'-wave amplitude, which is assumed to be 
the same for all solutions. 

The unambiguous determination of the two photon width of each resonance requires 
the corresponding amplitudes to be continued to the pole on the nearby unphysical sheet 
of the energy plane. This is done by following the procedure set out in Eqs. (42-49) 
in Ref. [16]. The residue of the relevant partial wave 77 — > w amplitude at the reso- 
nance pole determines the product of couplings g^g m (or equivalently the functions an j A of 
Eqs. (5,6) fix the ratio g~,/g n ). The residue of the corresponding amplitude, TJ^tity — > nn) 
given here in Eq. (5), then determines g\. From these the two photon width is given by 

where a ~ 1/137 is the QED fine structure constant. The resulting two photon widths 
for the a, /o(980) and ^(1270) are listed for Solutions A and B in Table 3. This is the 
primary output of this analysis. 

Our global fits determine the two photon width for the /o(980) to be 415 eV for the 
favoured solution A. However this lies in a range from 540 eV (for a solution qualitatively 
very close to A), down to 96 eV (as illustrated by Solution B). This is to be compared 
with various calculations dependent on the essential composition of this interesting state. 
A simple ss structure gives ~ 200 eV as found by Barnes [31]. His calculation [35j of a 
KK composition in the molecular model of Weinstein and Isgur [36.] gives ~ 600 eV. In 
contrast, more recently Hanhart et al. [37] predict 220 eV in a i^i^-molecular picture. 
This is comparable to the much earlier calculation of Achasov et al. [38J, who find, with 
a qqqq composition, a prediction of ~ 270 eV, using Achasov's kaon loop model. Despite 
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Figure 21: Argand plots for the 77 — > nn I = and 1 = 2 S- wave amplitudes. 
The normalisation of these amplitudes, Fjr° =0 , is defined by Eqs. (1-3). On 
the left are the 1 = S-w&ve for solutions A and B. The labels and bigger 
dots mark the energy every 0.1 GeV. The smaller dots are the intermediate 
energies every 25 MeV. The amplitudes A and B are very similar below 600 
MeV, by construction, cf. Fig. 19. As seen the amplitudes move particularly 
fast between 950 and 1000 MeV because of the /o(980), that in A showing a 
somewhat bigger loop. The amplitudes display the expected "kinks" at K + K~ 
and K°K thresholds. On the right is the common 1 = 2 S-wave with the 
dots labelling the energies every 0.1 GeV. 
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the claims of differing compositions, the essence of these latter two calculations is the 
same. If the / (980) width is indeed ~ 415 eV as in our favoured solution, this does 
not point unambiguously at any one particular modelling. However, as we have noted, 
solutions with smaller values of the 77 coupling are possible — at present at a lower level 
of statistical probability, though nevertheless acceptable. Thus the Achasov et al. [3H] 
and Hanhart et al. |37J prediction may yet be realised. Clarification requires precision 
77 — > 7r°7r° results, as we discuss in detail below. 

For the / 2 (1270) the helicity zero fraction is given — though determined from the 
pole residue, it accords very closely with the ratio of cross-sections shown in Fig. 20. 
Folded into the error quoted on the /2(1270) is the overall systematic uncertainty in each 
experiment. Though these are typically 10-12%, the tension between the different datasets 
and their systematic errors means that a typical fit results in a smaller uncertainty on 
any given solution. 

The two photon width of the a is obtained by analytically continuing the 7 = 5- 
wave amplitude, shown in Figs. 19, 20, to the <r-pole. The irn amplitudes we use in 
this analysis have this pole located at the position found by Caprini et al. [39J, viz. at 
E = 441 — z272 MeV. The analytic continuation is performed using the twice subtracted 
dispersion relation described in Ref. [23], but including not just it, p, cu-exchanges in 
the t and w-channels (Figs. 4, 5), but also the a\ and 61, which the work of Ref. [2JJ] 
indicates are not insignificant. This determines the ratio of the coupling a — > 77 to that 
for o — > irn, i.e. g-y/g n - The "radiative width" of Eq. (7) requires the o — > rnr coupling to 
be determined at the pole too. This value was estimated in Ref. |25j to be 0.55 GeV (see 
Eq. (4) of Ref. [23] for the definition), while Leutwyler has given a smaller preliminary 
value of (0.47 ± 0.02) GeV [27j. Using the latter value gives the results listed in Table 4, 
which are inevitably smaller than that found in Ref. [23] closer to those in Ref. [28J. 

In Fig. 22 are shown the total cross-sections for 77 — > tt + tt~ and 7r°7r° for both So- 
lutions A and B evaluated in 20 MeV bins. If the / (980) radiative width is "large" as 
in Solution A, with a width of 415 eV, then the peak in the 950 MeV region becomes 
even more pronounced in both the charged and neutral pion data if the full angular range 
were accessible. In contrast, a much smaller 1 = S- wave component is compatible with 
the Belle charged pion data over its limited angular range, if the peak near 950 MeV is 
produced by sizeable S — Do wave interference. Such interferences disappear over the full 
angular range when the partial wave components become orthogonal. This explains the 
difference in cross-sections between Solutions A and B in Fig. 20, where the contribution 
from individual partial waves with 1 = and 2 are shown. 

Precise measurement of the 77 — > 7r°7r° channel can distinguish between these two 
classes of solution. Not only are the total cross-sections different as displayed in Fig. 22, 
but importantly the angular distributions shown in Fig. 23 are quite distinct particularly 
at large z = cos 8*. With data promised from Belle in the near future, this ambiguity 
may well be resolved. 
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Figure 22: Total cross-sections for 77 — > 7r + 7r~ and 7r°7r° for Solutions A and B 
evaluated in 20 MeV energy bins with the bin centres connected by continuous 
lines. 
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Differential cross-section yy— 071° 71° 
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Figure 23: Differential cross-section for 77 
evaluated in 20 MeV energy bins centred on the energies specified in GeV. Even 
if the forthcoming precision results from Belle cover only | z = cos8* \ < 0.8, 
it should be possible to differentiate between these solutions. 
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Table 3: Two-photon partial widths in keV corresponding to the two classes 
of solutions as determined from their pole residues. For each the fraction of 
the width provided by helicity zero is given: for the scalar resonances, this is, 
of course, 100%. The error on the sigma's two photon width is largely due to 
the uncertainty in the a — > ixix coupling at the pole, discussed in Refs. [26J [28] 
and in the text here. 







Jo^oUUJ/a 




Pole positions (GeV) 


1.276 -^0.094 


0.441 - i0.272 


1.001 -z0.016 


r(/j — > 77): solution A 


3.14 ±0.20 


3.1 ±0.5 


0.42 


helicity zero fraction 


13% 


100% 


100% 


r(/j — > 77): solution B 


3.82 ±0.30 


2.4 ±0.4 


0.10 


helicity zero fraction 


26% 


100% 


100% 



5 Discussion 

An essential prerequisite of a meaningful Amplitude Analysis of the process 77 — > nn is 
differential cross-section data covering the widest possible angular range for both the 
7T + 7r - and 7r°7r° channels. While the Belle collaboration provides data of enormous statis- 
tics on the charged pion channel, this only covers 60% of the the range in cos#*. To 
make an Amplitude Analysis feasible, the data on the neutral pion channel from Crystal 
Barrel [TH [12], covering 70-80% of the angular range, with limited statistics (see Table 1, 
and Figs. 7, 11) are given sizeable weight, so that they play a role of equal importance to 
the several datasets on the 7r + 7r~ mode. The resulting solutions are of good quality, as 
indicated in Table 2 and Figs. 7-18, following all datasets adequately. While Solution A 
with a relatively strong /o(980) signal corresponding to a two photon width of 415 eV is 
preferred in terms of \ 2 (Table 2), partial wave solutions with a range of /o(980) com- 
ponents are possible. As illustrated by Solution B, a two photon width of the /o(980) 
as small as 100 eV is at the extreme of acceptability. This is made possible because the 
limited angular coverage allows the data to be described by a smaller 5-wave, but larger 
Dq wave, as shown in Fig. 20. At the other extreme a solution, with the qualitative fea- 
tures of Solution A but much higher x 2 , has a radiative width for the /o(980) of 540 eV. 
Precision results on 77 — > 7r°7r° should resolve this remaining ambiguity, as illustrated in 
Figs. 22, 23. 



31 



Achasov and Shestakov [40J have analysed the Belle integrated 7r + 7r~ cross-section in 
terms of model amplitudes in which key resonances have both a direct and meson loop 
contributions. While they find adequate fits to the charged pion data, their amplitudes 
fail to reproduce the 7r°7r° cross-section of Crystal Ball above 800 MeV. In contrast, our 
Amplitude Analysis treats both charged and neutral channels simultaneously and fits 
adequately not just the integrated cross-sections, but the angular distributions as well. 

In the simplified analysis of Ref. [40] . the Belle data lead to a somewhat larger 
/2(1275) — > 77 width than previously found. This is at first sight not surprising look- 
ing at the comparison of experimental data shown in Fig. 2. The PDG quotes [T7j 
T(/2(1275) — > 77) = 2.60 ± 0.24 keV largely based on the earlier (pre-Belle) Ampli- 
tude Analysis of Ref. [16], which gives 2.84 ± 0.35 keV. Our presently favoured solution 
A, in which the f 2(1275) is produced predominantly with helicity two, gives a compa- 
rable result of T(/2(1275) — > 77) = 3.14 ± 0.20 keV. In contrast, Solution B with its 
large D component, the ^2(1275) width rises to 3.82 ±0.30 keV. Once again forthcoming 
77 — > 7T 7r° data have the power to reduce these uncertainties dramatically. We eagerly 
await the completion of the analysis of these data. Their inclusion in an Amplitude 
Analysis will hopefully lead to a consistent set of two photon couplings for the low mass 
isoscalar states: /o(980) and /2(1270). Two photon couplings are a key window on the 
detailed structure of these low mass isoscalar states. Results from Belle make definitive 
conclusions for these within reach. To reduce still further the uncertainties in the 77 
width of the a requires precision charged and neutral pion data between threshold and 
700 MeV [3J. The introduction of appropriate taggers in an upgraded DA$NE machine 
at Frascati [4T] may well make this feasible too. We may then be certain which of the 
possible compositions [5J IH [5J these scalars have. 
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